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Abstract 

The non-divergent diagrams describing two-gluon exchange and 
annihilation between quarks and antiquarks are calculated in the 
Feynman gauge, based on quantum chromodynamics in a spherical 
cavity. Using the experimental N, A, f2, and p masses to fit the 
free parameters of the M.I.T. bag model, the predicted states agree 
very well with the observed low-lying hadrons. As expected, the 
two-gluon annihilation graphs lift the degeneracy of the tt and rj, 
while the p and uj remain degenerate. Diagonalizing the rj — rj 
subspace Hamiltonian yields a very good value for the mass of the 
r\ meson. 
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I. INTRODUCTION 



Cavity quantum chromodynamics, i.e. quantum chromodynamics with field 
operators obeying the linear boundary conditions of the M.I.T. bag model 
on a statics sphere (CQCD) [[IJ-0], is a consistent relativistic quantum field 
theory on its own. Indeed, CQCD may be expanded perturbatively M, and 
the diverging Feynman graphs have been shown to be renormalizable 
e.g. in the MS scheme. Taking into account the quadratic boundary condition 
of the M.I.T. bag model, the properties of the low-energy hadrons have been 
calculated successfully to order a s [0-0. 

However, this good agreement has been obtained by neglecting the self- 
energies of the quarks. More recently, the self-energies have been calculated, 
but unfortunately they turned out to be quite large thus spoiling the 

good agreement to order a s . One has therefore argued that the self-energies 
should perhaps be discarded, because the boundary conditions already account 
for at least part of them. 

Moreover, it has been pointed out that the perturbative expansion of CQCD 
in terms of a power series in the large strong coupling constant may be ill- 
defined. However, one can argue that the actual value of a s is immaterial, as 
long as the spectrum fits, order by order, the same hadronic states, N, A, Q, 
and p, thus fixing the free parameters of the model. Furthermore, it seems 
obvious that perturbative expansion in terms of cavity modes, rather than 
plane waves, is a much better starting point for the description of finite size 
hadrons. 

In spite of all these arguments, it is surprising to note that there was in 
the recent past little enthusiasm to develop CQCD beyond first order to check 
whether these expectations are actually fulfilled or not. It is, therefore, the pur- 
pose of this paper to clarify this issue, by calculating all non-divergent graphs 
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to order a"l in CQCD that can be written in the form of two-body operators, 
and to compare the resultant low-lying hadron spectrum with the experimental 
data. We are concentrating only on the two-body interactions in this paper. 
A calculation of the nondivergent three-body interactions for massless quarks 
has been performed showing that they are of much less importance than the 
two-body interactions to order a^. These interactions leading to a three-body 
force are not included in the following. 



II. SECOND ORDER ENERGY SHIFT 



Using the symmetric form of the Gell-Mann and Low theorem due to Sucher 
@, we can write the energy shift of an eigenstate \4> k ) of the unperturbed 
Hamiltonian as 

ied{^) k \S^ k ) c /drj 



E k -E? = lim 



x-i 2 



(1) 



{(j) k \S^\(j) k ) c 

where the subscript c indicates that only connected diagrams are included. The 
adiabatic S-matrix may be expanded as 



^ = i + E ^ (n) 



(2) 



n=l 



with 
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Here, the parameter e > introduces the adiabatic switching on of the inter- 
action Hamiltonian H[ nt (t) in the Dirac picture. Expanding AE up to second 
order in a s , we find that the only non-divergent terms contributing to AE are 
(S £ ^) c and (S e ^') c , and thus AE reads to that order 



AE = lim ie 



(S^) C + 2(S^) C 



(4) 
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The first term gives rise to the well-known one-gluon exchange and annihilation 
graphs m, while the second contains the two-gluon exchange and annihilation 
graphs which we would like to evaluate in this paper. 

Stoddart et al. |IJ have calculated the two-body operators in second order 
including six of the 24 possible time-orderings for both the two-gluon exchange 
and annihilation graphs. Here we take into account all time-orderings, using 



the results of ref. |10[ as a check. In fact, all the results of ref. |10[ can be 
reproduced with our method. 

As an example, let us calculate the 'straight' two-gluon exchange graph 
(Fig. |I](s2gx)). The energy shift due to this diagram is given by 



AE = lim 2ie / ax\ d 4 X2 0X3 d 4 Xi e 

X (fa 



-e(M + M + |t3|+M) 



T 



i9^4j>) f^yA^) (ig^ijj (ig^4 d $ 



/XI \ 

I I 



'X% 



X4-I 



L 



J 



(5) 

The time integrals are readily evaluated and the space integrals combine into 
quark-gluon vertex integrals defined by 



Qnn' = i d xu n (x)^^u nl (x)a^{x), 



(6) 



where u n {x) and a^^x) are the quark and gluon cavity modes, respectively 



111. The result is 
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Here, we have used the quark and gluon propagators in the Feynman gauge 



11]. The subscripts of the quark creation and annihilation operators stand for 



the color, flavor, and orbital quantum numbers, respectively. Repeated indices 
on the Gell-Mann matrices indicate a summation according to the Einstein 
convention. The quark-gluon vertex integrals are combined into the function 
Qs which depends on all the quark and gluon quantum numbers involved in 
the process, 

S - Wfq Hf'q'^qi ^q'i' ~ W fq Wf'-q'Wqi W-q'i'&ll 

rymYj nmE fym'YI fym'YI 771 , nmE nmE nrn'S' nm'E' 771 
~~ tyf-q Wf'q' ^4-qi ^£q'i' -^HI T tyf-q Wf'-q' V -qi W-q'i' -^IV- 

The terms E\ to £"iv, given in Appendix ^|, are sums over the energy denomi- 
nators that arise from the time integrations. 

The energy shift can be interpreted as a two-body operator V in second 
quantization sandwiched between the states \4>k)- This operator can be trans- 
lated into first quantization, yielding 

^12 = f "12(A), (9) 

where A stands for the total exchanged angular momentum between the quarks. 
F12 denotes the two-body operator F 12 = F\ ■ F 2 . 

Expanding the quark-gluon vertex integrals into angular and radial parts 
(S^f), we arrive after some algebra at 

ma) ^-E£ EE CLCT (-D j ' +Jj+a 



EE' jqj q l J1J2 A ^fi^-Rfi^z-R 
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Here, the summation over the radial quantum numbers of the intermediate 
quarks and gluons have been omitted for simplicity. The factors G and S are 
defined as 



Gj lJ2 (f, q, i) = JfJifi -h-h (-1)'/+*"* 
and 
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Since we are interested in the interactions of quarks and antiquarks 
ground-state, we can restrict ourselves to j n = | with n = i, i', f, /', for 
Eq. fliTf ) simplifies to 



(12) 
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The ^12 can now be determined numerically. We note that our method differs 



from that of ref. [ID in which the sum over all intermediate states coupled to 
exterior states with j = or j = 1 was calculated, while here we evaluate the 
coefficients of the two-body operators directly. 



III. RESULTS AND DISCUSSION 



The two-body operators for the energy shifts due to the interactions not 
discussed above are may be derived in the same manner. In the following, we 
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take all the two-body interactions into account. The two-body operators of the 
energy shifts can be written as 

V£ sx) = a s C[l gx) (A^ + B (1 ^ S 12 ) (15a) 

for the one-gluon exchange graph, 

= a 2 s Cf 2 gx) (A^ + B<**> S 12 ) (15b) 

for the straight and crossed two-gluon exchange graphs, and 

V g ga) = a 2 s (± - T 12 ) C[l sa) (A {2 ^ + S 12 ) (15c) 

for the straight and crossed two-gluon annihilation graphs. The two-body 
operators C\ 2 and the coefficients A and B are found in Table |. V\ 2 is given in 
natural units he/ R, R being the cavity radius. Si 2 , T 12 and F 12 stand for the 
product of the spin, isospin, and color operators of the two external particles, 
respectively. 

The one-gluon annihilation graph only contributes to the energy shifts of 
quark- ant iquark pairs with the quantum-numbers of a gluon. These do not 
occur in the hadrons considered here. However, for completeness we quote the 
resulting two-body operator as well, 

n ( 2 lga) = 0.187505 a s (± - T 12 )(F 12 + §)(S 12 + |). (16) 

At first sight, the results for the two-gluon exchange and annihilation seem 
to disagree with those previously obtained by Stoddart et al. [|10j, but since, 
in contrast to ref. [[H]], we have included all the time-ordered graphs, this is 
not surprising. If the coefficients for the two-gluon exchange diagrams are 
calculated to the same maximal energy including the time-ordered diagrams 
calculated in ref. ]l0| only, we are able to reproduce those results. For the 
two-gluon annihilation graphs, we obtain the same result for the energy-shift 
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of a quark-antiquark pair in a state j = 0, even though our coefficients differ 
from those of ref. for the unphysical case j '• = 1. 

Using the mass formula of the M.I.T. bag model |T|-§[, 

M = Uatt B) 1/4 n 3/ \ (17) 

o 

where B is the bag constant and Q is the total energy of the state in units of 
hc/R, we arrive at the mass spectrum for the low-lying hadrons (Fig. The 
total energy Q is of the general form 

n = Nuj - Z + V x a s + V 2 al (18) 

where N is the number of quarks, u the single particle energy, Z the vacuum 
energy, V\ and V2 are the first and second order energy shifts, respectively. 
Consistent with the philosophy mentioned above, we neglect possible contri- 
butions of order a s and al to uj and Z. The masses of the less problematic 
hadrons, i.e. the nucleon, the A-resonance, and the p-meson, are used to fix 
the parameters a B , B and Z. The fi^-hyperon fixes the mass of the strange 
quark. The parameters obtained are 

a s = 1.008, (19a) 

Z = 1.471, (19b) 

= 158.0 MeV (19c) 

m s = 1.445 fin" 1 = 285.1 MeV. (19d) 



The fact that including higher order graphs lowers the value of a s from 2.2 |]I2 
in first order to just over unity in second order is very encouraging. 

The two-gluon annihilation graphs lift the degeneracy of the ir and the r] 
mesons while keeping the p and uj degenerate. As we are able to calculate the 
diagonal and off-diagonal matrix elements due to two-gluon annihilation, we 
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may follow the procedure outlined by Donoghue and Gomm |13| and diagonal- 
ize the Hamilton matrix, which in the basis (uu + dd) / \/2 and ss reads 



n 



f fill fii 2 ^ 



^21 ^22 



(20) 



with 



fin = 2u u -Z + a s V^u ^u) + a 2 s {V 2je {u -> u) + 2V 2 , a (u -> u)) , (21a) 

fii 2 = V2a 2 s V 2 , a (u -> s) , (21b) 

fi 2i = v^a 2 K 2 ,a(s -> u) , (21c) 

fi 22 = 2cu s - Z + a s K(s -> s) + a s 2 (V 2 , e (s -> s) + V 2i „(s -> s)) . (21d) 

Here, u denotes a massless quark, s a strange quark, and the labels e and a 
refer to interactions due to gluon exchange and annihilation. We differ slightly 
from the procedure of Donoghue and Gomm JTS] , as they use the experimental 
values of the pion and kaon masses in their calculation. Deriving the masses 
of r) and rf in the framework of the M.I.T. bag model by diagonalizing the 
submatrix fi and inserting the eigenvalues into Eq. (jI7|), we obtain 



m„ = 534 MeV, 



m v > = 742 MeV, 



(22a) 
(22b) 



with a mixing angle of 33°. The resulting 77 mass is very close to the experi- 
mental value of m v = 547 MeV, while the rj' mass is too low compared with 
m v i = 958 MeV. However, here we have only taken into account contributions 
from massless and strange quarks to the mass of the rj'. The rf meson could 
contain a significant cc component which might well account for the remaining 
difference, since heavy quark masses have an important influence on the energy 
eigenvalues. 
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As one can see in Fig. the fit is very good for the masses of the baryons. 
The splitting between the A and the S is much better than in the calculation 



to first order in a s [|T^] . However, the predictions for the vector mesons </> and 
K* are still too high, while for the pseudoscalar mesons, tt, K, rj and rj', they 
are too low. 

The fit only takes into account the two-body diagrams calculated in this 
paper, but none of the divergent graphs to order a^, i.e. those that have a 
loop on any of the external or internal lines or vertices in Figs. |l](lgx) and 
[l|(lga). These divergent graphs should be included to make the second order 
energy shift gauge invariant. However, even though the cavity renormalization 
techniques have been developed , the numerical effort required to calculate 
such graphs is considerable^ and beyond the scope of this paper. In any case, 
these contributions are of the same form as the two-body operators for one- 
gluon exchange and annihilation. They will merely renormalize the first-order 
contribution by a factor. Thus the part which is orthogonal to the first-order 
contribution will still be gauge independent. The fit has been found without 
the introduction of this factor. 



IV. CONCLUSION 

We have calculated the non-diverging Feynman diagrams up to second order 
in a s in the framework of cavity quantum chromodynamics and obtained the 
two-body operators for the energy shifts. The free parameters of the model, the 
strong coupling constant a s , the zero-point energy Z , the bag pressure B, and 
the mass of the strange quark m s , are fixed by four out of the sixteen available 



The regularization procedure relies on delicate cancellations between oscillating 
terms. 
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masses of the hadrons containing only up, down, and strange quarks coupled 
to total angular momentum of up to 3/2. The calculated masses agree very 
well with the experimental data. It is especially noteworthy, that the mass of 
the ?7-meson can be calculated by diagonalizing the 77-7/ subspace Hamiltonian 
without introducing any further parameters into the model. 

The inclusion of second order interactions leads to a much smaller value of 
a s 1) than only considering first order terms (ct B > 2). Self-energies will be 
included in future work, but at the moment we neglect them as they might be 
taken care of by the boundary conditions already. 

There is still a lot of work ahead, to complete the calculations in cavity 
QCD to order a^. This includes the evaluation of divergent diagrams, but 
also the calculation of other observables than the masses, e.g. the magnetic 
moments, charge radii, and the ratio QaIqv- 
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APPENDIX A: "STRAIGHT" TWO-GLUON EXCHANGE 
ENERGY-DENOMINATORS 
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FIGURES 
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FIG. 1. Finite Feynman diagrams of first and second order in a s . 
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FIG. 2. The hadron spectrum. 
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TABLES 
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= m s 
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0.052685 
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-0.084483 


mi 
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-0.457013 


0.551227 


0.009522 


-0.017055 


0.036234 



TABLE I. The numerical results for the coefficients of the two-body operators 
as defined in Eqs. (|l^). The abbreviations stand for one-gluon exchange, straight 
two-gluon exchange, crossed two-gluon exchange, straight two-gluon annihilation, 
and crossed two-gluon annihilation, respectively. The second row gives the appro- 
priate color factors. 
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